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The parton momentum density distributions of the proton were obtained from a 
NLO QCD analysis of HERA and fixed target structure function data. The result- 
ing parton distribution set includes the full information on errors and correlations. 



1 Introduction 

Standard sets of parton densities i are widely used to calculate hard scat- 
tering cross sections in hadron-hadron and lepton-hadron collisions. However, 
none of these sets give the errors on the parton densities which tend to domi- 
nate the uncertainties on the predicted cross sections. 

To make a parton distribution set available which includes the full infor- 
mation on errors and correlations^ we have performed a NLO QCD analysis 
of HERA and fixed target structure function data. □ In this contribution we 
describe how the experimental errors were propagated in the QCD fit. 



2 Error propagation 

In the fit, the correlated experimental systematic errors were incorporated in 
the model prediction of the structure functions: 



F,( P , S ) = Fr(p)(i+E^ri (i) 



where F^ CT> {p) is the QCD prediction and A^ st is the relative systematic error 
on data point i stemming from source A. The fitted parameters p describe the 
parton densities at the input scale Qq and s denotes the set of systematic 
parameters. Assuming that the s\ are uncorrelated and Gaussian distributed 
with zero mean and unit variance, the x 2 was defined in the usual way and 
two Hessian matrices M and C were evaluated at the minimum \ 2 : 

1 d 2 X 2 1 d 2 X 2 

M X fj, = -- — q — , Cxf_, = - j — 5 — - (2) 

2 dpxdp^ * 2dp\ds fl 



"Available from ittp://www. nikhef.nl/user/h24/qcdnum 
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The statistical and systematic til covariance matrices were calculated from 

v stat = M -l^ v sy S t = M -1 CC T M -1 (3) 

The error on any function F of the parameters p is, to first order, given by 

( a ^) 2 =ee!^# (4) 

^ dp\ dp,j, 

where V is the statistical, the systematic, or if the total error is to be calculated, 
the sum of both covariance matrices. 

Three additional sources of error were considered in the analysis: (i) Er- 
rors due to the uncertainties on the input parameters (a s etc.). The error 
bands are defined as the envelope of the results from the central fit and two 
additional fits where each input parameter was lowered or raised by the error; 
(ii) An 'analysis' error band is defined as the envelope of the central fit and 10 
alternative fits where, for instance, the cuts on the data were varied; (iii) The 
renormalization and factorization scale uncertainties, obtained from fits where 
both scales were independently varied in the range Q 2 /2 < /.i 2 < 2Q 2 . 

In Fig. ^ we show the parton densities obtained from this analysis (left 
hand plot). The error bands correspond to the quadratic sum of all errors 
except the scale uncertainties. The relative error contributions to the gluon 
density and the singlet quark density are shown in the right hand plot of 
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Fig. [j]. It is seen that the analysis error band is small. For the gluon density 
the remaining contributions to the error are roughly equal in size whereas for 
the quarks it turns out that the scale uncertainty is the largest source of error. 



3 Parton distribution set 

Stored on a computer readable file are the statistical and systematic covariance 
matrices, the parton densities /j, the derivatives dfi/dpx (both from the central 
fit), the results from the systematic fits (where the input parameters or scales 
were changed) and the analysis error band. The kinematic range covered by 
the parton densities is 9 x. 10~ 4 < x < 1 and 1 < Q 2 < 9 x 10 4 GeV 2 . 

A computer programla gives fast access to these results and provides tools, 
which make use of Eq. (0) , to propagate the statistical and systematic errors 
to any function F of the parton densities. As an input to the program the 
user should provide a calculation of the derivatives dF/dp\. Let us take as an 
example a hadron-hadron cross section which can be written as a convolution 
of the parton densities and a hard scattering cross section, 



y^/i® fj®<rjj- (5) 



To calculate the error on a it is sufficient to provide a function which computes 
the derivatives 



da x - 



op\ dpx 



(6) 



This calculation is straight forward since the fi and the df jdpx are available 
from the input file. 

Finally, we remark that the errors from a QCD fit are not determined by 
the experimental errors alone but also depend, and maybe quite strongly, on 
the assumptions made in the analysis, in particular on the parameterizations 
chosen for the parton densities at the input scale Qg. 
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